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ABSTRACT 


The differential equations .of motion, and boundary conditions, describing 
the flap- lead/ lag-torsional motion of a flexible rotor blade with a pre-cone 
angle and a variable pitch angle 6(x,t), which incorporates a pre-twist 
are derived via Hamilton's principle. The equations are valid tor. both ex- 
tenslonal and inextenslonal Wlades, and the meaning of inextenslonallty is 
formally discussed. The equations are reduced to a set of three Integro partial 
differential equations by elimination of thv<! "extension"' variable. Both cases of 
hover and of forward flight are addressed. The generalized aerodynamic forces 
are modeled using Greenberg's extension of Theodorsen's strip theory. After the 
equations of motion are obtained, they are systematically expanded into poly- 
nomial non-llnearltles with the objective of retaining all terms up to third-degree 

so that the influence of such terms on the motion of the system may be evaluated. 
The blade Is modelled as a long, slender, initially straight beam of Isotropic 

Hookean material. Offsets from the blade's elastic axis through its shear center 
and the axes for the mass, area and aerodynamic centers, as well as radial non- 
unlformltles of the blade's stiffnesses and cross section properties are taken 
into account. The effect of warp of the cross section Is also Included in the 
formulation. 
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CHAPTER I - INTRODUCTION 


An Importanc problem assoclaced with helicopter dynamics Is Che decermlns' 
cion of Che dynamic response and aeroelasclc scablllcy assoclaced with Che rotor 
blades. Considerable attention has been given to rotary-wing aeroelastlclty pro-’ 
blems, as evidenced by Che numerous publications In >nls area, e.g. (1-60]. As 
a rejult of this effort, It Is now recognized that such problems arc Inherently 
non-linear and thus special attention should also be given to the development of 
consistent mathematical models to address them. Due to their flexibility, rotor 
blades are able to bend in two principal directions as well as twist. This fact 
alone introduces several important non-linear terms In the differential equations 
of motion. It Is also well known that considerable complexity Is further intro- 
duced due to the effect of rotation, through the addition of centrifugal. Inertia 
and Coriolis forces, and of the aerodynamic forces. Comprehensive reviews of the 
subject have been written by several Investigators, e.g., [1-3]. The work reported 
In [1] by Loewy reviews the literature up to 1969, while the recent review presented 
by Friedmann (3) Includes most of the work done up to 1976/1977. Several aspects 
of helicopter modelling have been discussed by Ormiston (4). 

In 1938, Houbolt and Brooks (14] derived a comprehensive set of linear 
differential equations to describe the coupled f lap-lead/lag-torslonal (or. In 
short, flap-lag- torsional) dynamics of non-uniform rotating pre-twisted blades. 
Expressions for the aerodynamic loads were not specified in [14]. Further work by 
several Investigators established that non-linear terms not considered in [14] play 
a fundamental role in the response and the stability of both hingeless (cantilevered) 
and articulated (pinned) rotors. The non-linear equations of motion are quite 
complex, and efforts are still devoted today toward their formulation. Hodges and 
Ormiston [15] analyzed the coupled flap-lag behavior, which can exhibit insta- 
bilities, by restricting themselves to torsionally rigid blades. A set of Integro 
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partial differential equations of motion with quadratic; non-linearities were 
developed taking into account the effects of the blade’s pitch angle and of a 
small pre-cone angle. A study of the stability of the flap-lag motions rf a 
hlngelesa rotor In hover was then performed by making use of such equations. A 
comparison between the results of the stability analysis for elastic and rigid 
blade models wae also presented In [IS]. Work with rigid blade models has also 
been done by several other investigators, e.g. (22-31)* 

Hodges and Dowell [32] developed a comprehensive set of differential equations 
with quadratic non-llnearltles for the flap-lag-torslonal dynamics of rotating non- 
uniform extenslonal blades. The presence of pre-twist and of a small pre-cone angle, 
and the effects of offsets between the mass and area centroids from the elastic 

axis through the shear center were considered. Three-axes Euler angles were used to 
describe the orientation of the principal axes of the blade's cross section — an 

approach that facilitates the development of the equations. This approach was also 
used by Peters and Ormlston [33] when considering the effects of the second-order 
non-linearities on the angle of attack of hingeless rotor blades. The equations In 
[32] were used in [34] to investigate the stability of a uniform rotor blade in the 
hovering flight condition. The equations of motion developed in [32] were also used 
by Dowell, Traybar and Hodges [36] to correlate the experimental and theoretical 
results obtained for the static deflections and bending natural frequencies of a 
cantilever beam with a tip mass, after the appropriate terms due to the Inclusion 
of the tip mass were added to such equations. Reasonably good agreement was 
reported for the cases where the basic assumption In [32] of retaining only quadratic 
non-linear terms in the equations could be justified. A similar study was also 
recently presented by Rosen and Friedmann [37]. 

The equations developed in [32,34] were further extended by Hodges [38] to 
include not only the effects of twist and pre-cone but also of droop and sweep 
angles, torque and blade root offsets, and pitch-link flexibility simulated by a 
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torsional root spring. 

The coupled f Isf^-lsg-pltch rssi^onss of hlngaless bladss for both cues of 
hover and forward flight were conaldared by Friedmann and Tong (40,41). A set 
of equations to describe the flap-lag-torslonal motions of a pre^cwlsted 
cantilevered blade In hover, considering root torsional flSKlblllty, waa aleo 
derived by Friedmann (42). The derivation presented In (42) follows Friedmann 

e 

and Tong's earlier work (40,41) and a modification of Houbolt and Brooke' work 
(14) to Include quadratic non-llnear terms. An ordering schema, also used by 
other Investigators, was employed to neglect several non'-llnear terms In a 
systematic way and obtain the equations with quadratic non-llnearltles. A revised 
set of non-llnear differential equations of motion was later derived by Rosen 
and Friedmann (48). The equations developed In (48) for pre-twlsted blades with 
a small pre-cone angle Include the effects of offset between the blade's area 
and mass centroids, shear and aerodynamic centers. As In (32), three axes 
Euler angles are used to describe the blade's orientation In space. 

Some of the work on the flap-lag dynamics of rigid hinged blades and hingeless 
elastic blades for both cases of hover and forward flight has been examined 
by Kara and Kvatemlk (50). For the rigid blade case, two types of hinges 
were considered, one for which the lead/lag hinge flaps with the blade, and the 
other for which the flap hinge also executes a lead/lag motion. Stability 
boundaries in the flap- lead/lag frequencies space were obtained. Indicating that 
such boundaries depend on the physical hinge sequence. This work was also 
reported in (30). The same authors also examined in (51) the flap-lag-torslonal 
kinematics of pre-twlsted blades using lag-flap-pitch and flap-lag-pitch three- 
axis Euler angle sequences to describe the orientation of the blade's principal 
axes. Here, however, a source of disagreement with the work of (32,33) arises 
when the third Euler angle — a final rotation about the blade's elastic axis 
— is equated In (51) to the twist angle of the blade rather than determining 
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Che blade's cvlaC from Che costponenC along che elaaClc axia of Che curvacure 
veccor (which la obcaloed as a funcclon of che pre-cwlsc, Che chree Euler angles 
and their spatial derivacivea) , as dona in (32,331. The kinematic expressions 
obtained in [51] are uaed in [52] to develop a set of non-linear differencial 
equations of notion with quadratic non-linearities for the forward flight 
condition using cvo-dimensional quasi-ateady blade element theory. The differ- 
ential equations for the two Euler angle sequence mentioned above are shown 
explicitly. Although not mentioned in [52], these equations, which represent 
che same physical problem, should be equivalent to each ocher. In [50,52] the 
authors make explicit use, in their derivation, of the blade's kinematic axial 
deflection, or foreshortening. 

The source of disagreements previously mentioned in regard to the formulation 
of che f lap-Iag>corslonal dynamics of rotating blades was recently clarified by 
Hodges, Ormlston and Peters [55]. They considered the kinematics of rotating 
extenslonal beams and developed expressions for che curvature vector and che Euler- 
angle transformation matrix for all possible six sequences of the three-axes Euler 
angles. By oK^ans of an appropriate change of variable, they showed that these 
sequences are equivalent to each other. In addition, they showed Chat the angle 
of twist, obtained by integration of Che twisting curvature^ does contain a higher- 
order integral term, and that such an angle is a quasi-coordinate much like che 
axial deflection used in [50,52] with Che foreshortening considered a priori. 

This angle of twist, first obtained by Peters and Ormiston in [33], is discussed 
by Nordgren in [61] and was independently obtained by Crespo da Silva and Glynn 
[62-71] in connection with the nonlinear flexural- flexural-torsional motions of 
inextensional beams. 

The non-linear differencial equations of motion for inextensional beams were 
derived in [62,68] by making use of Kirchhoff's kinetic analogy [72], three-axis 
Euler angles to describe Che orientation of che beam's cross sectional principal 
axes, ani. Hamilton's principle with che Lagrangian adjoined to the inextensional 
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constraint by a Lagrangian multiplier. Physically, the Lagrange multiplier repre- 
sents a force along the beam's neutral axis necessary to maintain Incxtenslonallty. 
This f /rce Is dependent on the deflections In a non-linear manner. The equations 
derived In (62,68] were then expanded Into cubic non-llnearltles and converted to 
a set of Integro partial differential equtitlons '63,68j. These were subsequently 
used to Investigate the non-linear response, and Its stability, exhibited by the 
beam under several conditions (63-67, 69-71] Including support asymmetry (64). 

The systematic formulation of (62, 63, 68] may be extended to address the rotor 
blade dynamics problem by Incorporation of effects such as blade rotation, pre- 
twist, pre-cone angle, and aerodynamic forces. 

An investigation of the response, and Its stability, of a dynamical system 
described by a set of nonlinear differential equations with quadratic and cubic 
nonlinearicies acting simultaneously has been presented by Nayfeh and Kamel, and 
by Nayfeh (73,74]. Although the problem considered in (73,74} Is also a complex 
one, the rotor blade aeroelastlc problem Is substantially more Involved. To 
this author's knowledge, a systematic Investigation of the Influence of the 
extensional and inextensional assumptions, and of the next order — cubic — non- 
linearities on the rotor blade aeroelastlc response is missing in the literature. 

In this report the work presented In [62,68,32] is extended with the objective 
of deriving the differential equations of motion for both extensional and inexten- 
sional rotor blades with pre-cone and a pitch angle 0(x,t) — which Incorporates a 

3 

pre-twist — including third-order (e ) non-linearities. The pitch angle 

8 may be included in several ways. One way, for example, is to start with an un- 
pitched blade and to impose the known pitch angle 8 about its longitudinal axis 
to obtain the directions of the principal axes of the undeformed but "pre-pltched" 
blade. Then a sequence of Euler angles can be used to bring these axes to their 
new orientation obtained after the blade's elastic deformations. Atother way, 
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introduced in [32,33], is to start with a set n^, ;») of non-principal axes 

of the pitched blade with, say,?^ along the blade's elastic, axis, and Cj) being 

two non-principal cress sectional axes at an angle -9 with the "pre-pitched" (but 
undeformed) blade's principal axes along its cross section. A sequence of three 
Euler angles is then employed to bring to their corresponding orienta- 

tion after the blade has undergone the elastic deformations. Since the blade's 
cross section is now, after the elastic deformations, in its "final" orientation 
in space, a further simple rotation 9(x, t) about the "final" position of the 

axis will then bring the triad into coincidence with the deformed 

blade's principal axes (^,n,?). With this approach, the pitch angle 0(x,t) appears 
in the differential equations simply as an aoditive term to the third Euler angle, 
as long as this angle is a rotation about the blade's elastic axis (32,33,55). Be- 
cause of its inherent simplicity and compactness, this is the approach that will be 
used here. 
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CHAPTER II - THE DYKAMIC SYSTEM AND BASIC ASSUMPTIONS 


Consider an Initially straight » but pre-twisted, rotot blade of length R 
when not elastically deformed, mass m per unit length, and of closed cross section. 
Its maximum cross section dimension is assumed to be much s'jaller than its length 
so that the blade may be approximated as a beam. A blade segment, both in its 
undeformed and elastically deformed states, is shown in Pig. 2.1. The (C* H* C) 

A A A A 

axes — with unit vectors (5,r|fC*Cxfl) — are the principal axes of the 

* 

cross section at the shear center C of the deformed blade; it is assumed that 

e 

the cross section is symmetric about the r) axis. The C axis is tangent at all 
times to the elastic axis of the blade. 

A A A A 

The (x,y.sz) axes, with unit vectors ( x,y,z-x x y ) shown in Fig. 2.1, are 
a set of rotating reference axes., the x axis is coincident with the elastic axis 
of the undefonned blade. These axes rotate in space with constant angular veloc- 
ity n about a direction perpendicular to the rotor hub. This is shown in detail 
in Fig. 2.2 together with their spatial orientation relative to a set of iner- 
tial (X, Y, Z) directions. As seen in Fig. 2.2, the orientation of these axes 
may be described by first aligning (x, y, z) with (X, Y, Z) and then performing 
two successive rotations. The first rotation T ■ fit ( where t is time) about Z 
brings the (x, y, z) triad to its new orientation (X,, Y, , Z,« Z), while a second 
rotation 3 about Che negative Y^ direction (i.e., a clockwise rotation) brings 
(Xi* Y^, Z^) to its final orientation (x, y, z). The angle 3 . which is the angle 
the undeformed blade x axis makes with the plane of rotation, XY, is the blade's 
pre-cone angle. This angle is taken to be a constant. The "blade root offset" e^ 
shown in Fig. 2.2 will be taken to be zero for simplicity. 

A 

The principal axes (n ) of the blade’s cross section, centered at C as 

e 

shown ia Fig. 2.1, make an angle 6(x, x ) with two cross section non-principal axes 
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also shown In that figure. It Is assumed that the angle 6(x,t) — the 
geometric pitch angle of the blade — is given as 


0{x,t) - 


c pt 


•2 

i - 1 




cos i T + 


is 


sin i T 


( 2 . 1 ) 


where 0^ is the collective pitch angle — a constant — , • pre-twist 

angle that may be incorporated to the blade, and 6^^ and 6^^, i ■ l,2,...,n, are the 
harmonic pitch components that may be introduced by a control system. When the 
blade is elastically undeformed the non-principal and cross section axes 
are parallel to the rotating y and z axes, respectively. 

During the elastic deformations, point C — the elastic center of the 
undeformed (but pitched) blade's cross section at location x«x shown in ?ig. 2.1 
— moves from its location (x*x, y^O, z*0) to Its new location whose 
coordinates relative to the (x,y,z) rotating a;^j» are written respectively, as 
Rx + Ru(x,t), Rv(x,t) and Rw(x,t). Here u, v and w are the (x, y, components, 

non-dimenslonalized by the undeformed blade's length R, of the elastic displace- 

\ 

ment vector of the blade's cross ^section elastic center. Clearly, 0 £ x £ 1 is a 
non-dimensional quantity. \ 

\ 

In genera^, each cross section of tn-e blade experiences the elastic displac'J- 

/ A 

'dents Ru(x,T'), Rv(x,x) and Rw(x,T) of its e'l|i8tlc center C and a rotation about 

* ^ \ ® * 

C^. The orienls^tlon of the principal axes (?,TtxC), through C^, of the blade's 






cross section at C • may be described by three successive Euler angle rotations. 

® . \ 

A set of three-axes Euler angles is used here for thisX^urpcse. We begin this 

^ A A AAA \ 

process by aligning (C, n» C) with (x, y, z) and then pern^rming the three 
successive rotations shown in Fig. 2.3 . The first rotation wibout z brings 


\ 


\ 


*Due CO warping of the cross section, the coordinate system (x, n, O Is in 
reality a non-orthogonal system if 36/ 3x ^ 0 since, in this case, C is not 
coincident with x. However, this effect Is quite^small [75, 76] and, for 
tnis reason, it will be neglected here. 
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A A A ^ . 

that unit vector triad to (Cj^t ■ *)• Th« aocond rotation 0 y(x,T) about 

A 

the negative direction (a clockwise rotation) of the new position of the 

A A A A ^ A 

n unit vector brings (Cj^. to (^2* " ’^ 1 * ^ 2 ^* ^ third rotation 

A A A A /\ 

9^(x,t) about ^2 * ^ brings this triad to the orientation (C* C^)* 

Thus, the sequence (9,6,6) takes a non'-prlncipal triad from Its Initial 

z y X 

orientation associated with the undeformed blade, and aligned with the rotating 
(x,y,z) reference axes, to Its *'final"orlentatlon associated with the elastically 

A A A 

deformed blade. As Indicated In Flg.2>l,an additional rotation of the (Ct H* C) 

A A A 

triad, which Is still aligned with (C« C 3 ) after the third Euler angle 

A 

rotation 6 ^(x,t), by an amount equal to the pitch angle 6 (x,t) about C will bring 
the blade's cross section principal axes to their "flnai* orientation In space. 

A A A A A /> 

The transformation matrix [T] between (x,y,z) and (^, Hi C) t defined as. 


m m 

A 



C 


X 

A 

n 

[T1 

A 

y 

A 


A 



t 


may be readily obtalnrid with the aid of Fig. 2.3. As found In [32, 33, 551 Its 


components t^^^ are given ns 



C0 

y 

c0 

z 




c0 

S0 







y 

z 

[T] - 

-c0, 

s0 

- 30, 

90 

c0 

c0. 

C0 - 80., 80 80 

1 

z 

1 

y 

z 

1 

z 1 y z 


30, 

s0 

- C0, 

S0 

c6 

-s0. 

C0 - C0- S0 80 


1 

z 

1 

y 

z 

X 

z 1 y z 



(2.3) 


where s9j^ and c9j^ denote, respectively, slnSj^ and cos0,^ (k-l,y,z), and 0^ Is given 
as 
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1 


0(X,T) 


(2.4) 


The advantage of Incorporating the pitch angle 9(x,t) in the way that is done 
here is that 9(x,t) simply shows up in the equations as an additive term to the 
third Euler angle 9 (x,t), regardless of the fact that the blade is given a 


known "pre-twist" general known pitch angle 9 (x,t). This approach 


was first introduced in [33] for pre-twisted blades. 


With the position vector of the elastic center C , relative to the hub 


center 0, given as 


r^“R [(x + u) x + vy+w^ 


(2.5) 


it follows that 


A 

^ e 

wcx 


■ (x-fu) x-fv y4^w z 


( 2 . 6 ) 


where 


( )‘^ A 3 ( )/3r . 


Equations (2.2) anc^ (2.6) imply that 


- c9^ c0^ - (x + u)^ -(!-♦- u^)x"*" 


(2.7a) 


t.^ ■ c9 90 ■ V • V X 

12 y 2 


(2.7b) 


* The symbol A is used here to denote ’’equal by definition”. The scalar 

quantity r denotes arc length, non-dimensionalized by the blade’s undeformed 
length R, measured along the elastic axis. 
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where ( )’ • 3( )/3x. Equations (2.7a-c) may combined to yield the 

^ t 9 ^ 

following expression for x In ccnns of u, v and w , 

* [(1 + u*)^ + v*^ + w*^l”^^^ 

Equation (2.8) will be used In the sequel to relate the partial derivatives 
( )^ to ( ) ' and u’ , v' , w* . 

By letting dots denote partial derivative with respect to Che non- 
dimensional time T "Qt, i.e. ( )* ■ 3( )/3T, the angular velocity 5(x,t) 

of the principal axis system with respect to the Inertial frame 

(X, Y, Z) Is obtained directly from Figs. 2.1, 2.2 and 2.3 as 

iU • Q[§ z-5 z 

z y 1 1 

- n [0, + (0 + c6)s0 + s0 c0 c0 ]C 

1 z y y 2 

+ ((0 + C0)80.C0 - 0 C0, - (c0, S0 + S0- S0 C0 )s6]n 

‘ z l y yl lz lyz 

+ n [(0 + c8)c0, c0 + 0„ s0, + (s0, 80, - c0, 80„ c0,)sB]; 

z iyyi i.z iyz 

AAA 

A [u)^ C + n + u)^ Cl 



By making use of Klrchhoff's kinetic analogy (72, Chap. 19), the components 


T 


Cp , C and C_ of the elastic axis curvature vector C(x,T) can be readily obtained 

^ n S 

from equation (2.9) as 


Rc - (e^ 2 - + e| t) 


- [e| + s6 it 19'*' s9, c0 - e”*" c0Jn 
‘ 1 z z 1 y y 1 


+ (9^ c0j^ c0y 0y s0j^]C A C n + tl 


( 2 . 10 ) 


Although a total of six dependent variables (u, v, w) and (0^, 0^, 0^) have 
been introduced, only four are needed. As seen from equations (2.7a-c), the 
angles 6^ and 6^ are related to the spatial derivatives u, v and w as 


tan 9 ■ v' /(I u* ) 

z 


(2.11a) 


sin 9y * w’x^ - w'[(l + u')^ + v'^ + 


(2.11b) 


The differential equations of motion to be developed here will involve 
the elastic displacements u(x,t), v(x,T) and w(x,T), and the third Euler angle 
9^(x,t) . To obtain such equations, the extended form of Hamilton's 
principle [77,78] will be used. For this, expressions for the variations of 
the kinetic and strain energies of the motion are needed. Expressions for these 
energies are developed next. 
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CHAPTER III - KINETIC AND STRAIN ENERGIES OF THE MOTION 


The poslclon of an arbitrary point P on the cross section through C , which 

6 

Is located by the vector r given by equation (2.5), Is shown In Fig. 3.1. 

Due to warping, P experiences a small axial displacement given approximately as 

+ ■' 2 
- (Rij;)(p^x )C, where O 1* the warp function — normalized by R — obtained 

by solving Laplace's equation for the cross section [79, 80]. In terms of non- 

dimensional coordinates n and C (along the r) and C axes, respectively), the 

position vector of P relative to the hub center 0 is then written as 

rp - r^ 4- R(n n + ; <;) - C 

- R((x + u + nt 2 j^ + 5t^j^ - ij)p^x'*’tj^j^)x 

+ (v + rjt 22 Ctj 2 - + (w + r)t 22 * ^^33 ’ ’^P^x'^'t^j) z] (3.1) 

The velocity of P relative to 0 Is then 

Vp - d?p/dt • Q((drp/dT)^^^^g + Z x tp) (3.2a 

where (dr_/dT)-^ ^ ^ is the velocity of P as seen by an observer fixed to the rotat- 
P Xf y p z 

AAA 

Ing reference frame (x,y,z). 

Equation (3.2a) yields 


16 



Vp/(nR) 


•fjv ^(x * u)cS - w«6 > ^^*22 ^ tj^cS - tjjxB) 

+Z(ij2 ’► t3iC0 - i33»S) - Cj^j^cB) - ly 


4'fw ■♦• v»B n(t 23 ■♦■ 


t22«B) 


•*■ C(t 


33 


•f tj2»8)*'<'p5* 


^^ 2 *® * 13 ^* 1 * 


(3.2b 


Assuming chst the vsloclcy of chs hub csncsr 0 Is constant , vhich Is trus 
for both cases of hover and of non-accelerated forvard flight, the kinetic 
energy of the motion of the blade is then 


T - 



// P V • V dn dc dx 
A 


(3.3) 


where o is the blade's material density at point P and A Is the undeformed blade's 

2 

cross section area, normalized by R . The mass per unit length of the blade, m, is 
defined as the following area integral 


m A 

a 


// P dn d; 
A 


(3.4a) 


and will be assumed to be a constant. This implies that P la not a function of 
X, that is, p*p(n,C). It will also be assumed here that the blade's cross 
section has material symmetry about the n axis, so that 


f/p C dn dj; - // p n C dn d; - 0 (3.4b) 

A A 
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f 




and that the warp function aatiafiaa the following relatione 


//pi|/dn dc - // pi|/ n dn dc - 0 

A A 


(3.4c) 


I With the aforementioned aaaunptiona , and defining the following croes oection 

Integrals 


m a ^ // pn dn dC 

A 


(3.5e) 


mj A // p dn dc ; mi ^ R^ // P dn d; 
'a ^ a 


(3.5b) 


^ ^n ^ 


(3.5c) 


the kinetic energy given by equation (3.3) may be written, after making use 


of the expressions for the elements tj^^ of the transformation matrix (T) given 


by equation (2.3), as 


T/(mfl^R^) " 1/2 / {(u-vc0)'^ 4> (v +(x, + u)c0 - ws0)^ + (w 4- vs0)^}dx 
0 


-n/2 / (.1^ ^n ‘^n 


0 f * * • • 

4- / e{(u-vc0) t22^3) + (w -4 vs6) 


23 "22^ 


(v + (x 4- u)c3 - WS0] (t22 + ^21*^^ " t22S0)}dx 


4- T*/(mn^R^) 


(3.6) 
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The terms not shown explicitly In equation (3.6) ere included in T . The 
expression for T is given In the Appendix. The quantities e and (a • Ct H» O 
defined by equations C3.5a-*c) are the blade's mass center offset from the elastic 
axis, normalized by R, and its distributed mass moments of inertia, normalized 
by mR^. 

To obtain the expression for the strain energy of the deformed blade, a 


strain tensor, represented here by a matrix [e] with elements is needed. 

As a measure of thw deformation of the blade, the equare of the distancee 
between two Infinitesimally close points on the blade, before and after the 
deformation, is used [80,81]. The nonlinear strain-displacement relations for 
the deformed blade are sought for the case of large displacements, and for this 
the strain components are formulated in terms of the increments dr, dr) and 
di;; for the deformed configuration. The strain components are still assumed to 
be small enough, however, so that Hooke's linear law relating the stresses 
at any point P of the deformed blade's material to the strains at P is still 
applicable. The matrix (£] of strain components is thus determined as 


dtp ‘ dtp - drpo 


drpQ - 2R^(dr dn d;] (el 


dr' 

dn 

d; 


(3.7a) 


where, as obtained from equations (3.1) and {' 1 . 10 ), 


PO 


u“v»*w"0 *0 

x 


R((x - «j^Q0’)x + (nQC0 - CgS0)y 


+ (hQsQ CqC0)z) 


(3.7b) 
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In equation (3.7b), Hq and denote, reapectlvaly, <(>, n and C when u • v • w • 

6^ ■ 0. Aa dlacuaocd in (32), the difference between the quantitiee and a 
(a ■ for the blade's croaa aection ia very eoall. In-plane croaa aeetion 

distortion will be neglected here, and tho approximations n, ; and ^ 

are adopted. 

From equations (3.1) and (3.7b) it follows that 
d?pQ/R • ((x*^- W- x”^)dr - 0'((3i(>/3n)dn 0tj^/30dCl}K 
+f-0'(na0 Cc6)dr (c0)dn - (a9)dOy 

+{0'(nc9 - Cs0)dr (S0)dn + (c0)dC)r (3*8) 

dtp/R • ((x -f u)^ X 4-v^y w"^* - ij;(p^x^)^<;]dr 

•Kn - C p^x'*’ 3ii;/3n)dn •«•(;-? P^x^' 3ij;/3c)d; 

+ x“’’vp^5: -f P^n + p^;) X (ri n *► C ; - ii»p^ x”^)dr (3.9a) 


Making use of equations (2.1), (2.2) and (2.7a-c), and defining, for the sake 
of compacteness, Xj^ - x, x^ • y and Xj ■ *, equation (3.9a) may be written, 
after performing the cross product Indicated and noticing that (p^x^)^ • x^(p^x^) ' , 
as 
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dFp/R - {(t^^ ^ (;p^ - np^)tj/ - 

■ ‘li 3<l»/9n)dn (tj^ - p^x'*’ 3tj^/aOd;)x^ 

Equations (3.7a), (3.8) and (3.9b) laBwdlataly ylald tha following axpvasalons 
for Che strain cooponants 

•(!•♦■ (CP^ - nP^)x'*’ - i|/x'*’(p^x'*’) + [(?*♦■ 'I'P^x"*’)^ (n “♦■ 'i'Pj^x'*’) ^ ) (p^x'*’) ^ 

- ((1 - ij^0")2 -*• (n^ + C^)9'^)x’*"^ 

2Ci 2 • - (; + ail'/sn) (p^ - e’)*'*' + ((nf^^ - ;p^)34'/3n - ij>pj.lp^x'^^ 

+ (p^x'*’(p^x'*’) - 0 6 Jx’^ijj 3ii//9n 

2 Ci 3 - (n - 91^/30 (P^ “ 9')x‘*‘ +((nP^ - CP^)9'I'/9C + <PP„lP^x'*"^ 

+ (p^x^(p^x"^) * - 9 9 ]x% 94*/ 

2 ej 2 • (9^/9n)^ 

2c 23 - ((P^x"^)^ - e'^] (9ii;/9n)(9ij^/90 

2 E 33 ■ [(p^x"*")^ - 9 (9i|;/9C)^ 


(3.9b) 


(3.10a> 

(3.10b) 

(3.10c) 

(3.10d 

(3.10e) 

(3.10f^ 
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Th« •Crain cnargy, U, la given by 


U • 


f'/// 

(blada) 


to 


11 U 


+ o 


22 


£22+033 C33 + 2(0 


12 '■12 


•f a 


13 *"13 


■f a 
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e23)l**n 


dC dx 

(3.11) 


where 0^^ are the screaies acclng on Che blade. For a linear and Isotropic 
elastic material, the stresses and the components of the strain tensor 
at any point P of the blade are relate 1 by Hooke *^s law as 


^11 " ^*^11 ‘ ''^‘^22 ^33^ ^^®y (3* 12a) 

?22 ■ {<^22 " ''’^‘^11 *^33^^'^y (3.12b) 

' 3 J * '”33 ■ “<‘’11 *'’22”/', 


£i 2 - (1 + V) 0 ^ 2 /^y " (3.12d) 

- (1 + V) Oj^j/Ey - o^2/^2G^) (3.12e) 


ej 3 • (1 + V) 023 /Ey - 023 /( 20 ^) (3.12f) 

where E^(x), v and Cg(x) • Ey/ (2(1 + v) ] are, respectively, the Young modulus, 

Poisson ratio and Che shear modulut; for the blade's material. 

With the assumptloii Chat the blade's maximum cross section dimension is much 
smaller than its length, the blade is approximated as a beam, and the normal 
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• tressea O22 ^33 Under this eseumptlon, equetlons ( 3 . 12 a-c) 

reduce Co 


"a = “ a '^ 


■22 


£33 5 - 


(3.13) 


and the strain energy U, calculated from equations (3«ll)v (3.12d*f) and (3.13), 
becomes 


“ ■ f' /' IK «(«, - 

x-0 ^ ^ 


np ^) x *^ +(? p ^ - np ^)^ x ‘’'^ + ;^) ( p | 


e ’ 2 ) K ^2 


+ f 


^( p ^ x "*") 


' ^ (P^ ^ pI)pI 


+(CP^ + np ^) p ^ x '*’^ - ( cp ^ - np ^) x '*'^ ( p ^ x ‘*')' J ^ 


x"*"^ - + 2tj^|e"x'*’^ - (p^x"^) x'*’ 




dn d^ dx 


* s ’ /’ IK 

^ A ® 


(; + 9 t |^/ 3 n )^ + (n - 


x»0 


+ 2p^(p^ 


|”(p^ - 0 ’)^ |^( 

- - CPp)94'/3C + tj^P^; 


r^2^ ■* 

O X 


+3 


-(c + 3>l^/3n) [jnp^ - Cp^)3i|;/3n - 


> dn dc dx + U 


(3.1A) 


where the term U is given in the Appendix, 
s 
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Aaaumlng chat the blade cross 

section Is symmetric about the n 

axis • and 


: the warping function 4^ is antisyvnnetrlc , it follows that 



// ; dn dc - //cn dn d; - 

A A 

//cn^ dn dc - //c^ dn dc - 0 
A A 


(3.15a) 

ff<p Cn dc - dn dC - 

A A 

//<l>(n^ c^)dn dc - 0 

A 


(3.15b) 

Non-dlmenslonalizlng the Young and the shear moduli as 

1 . 



E A Ey/(mfi^) ; C, A G^/(mQ^) 


(3.15c) 

defining the following non-dimensional cross section integrals, 

» 


Dr A G //[(c 3t|//8n)^ + (n - 

3ii</3c>^]dn d; 


(3.16a) 

^ “ A 




DAE //c^ dn d(: ; 

' A 

D A E//n^ dn dC 
^ “ A 


(3.16b) 

I A // (n^ + c'^)dn d; ; 

^ A 

A A //dn dc ; Ae.A 

A 

//n dn dc 

A 

(3.16c) 

B, A // (n^ + dn dC ; 

A 

B- A // n(n^ + c^)dn dc 
^ A 


(3.16d) 

^ dn dc ; 

A 

B 4 A B 2 - B 3 


(3.16e) 

A // ip^ dn d? ; 

A 

^ //? 'P dc 

A 


(3.16f) 

■Jp ^ (4> - i: 3ip/dO + (c + 3t/^/3n)c 3<|//3n]dn dc 

^ * A 


(3.16g) 
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(3.16h) 


Jr ^ // 34'/35 + (; + 3'i»/3n) ('!< - n 3<l'/3n)l dn dc 


the strain energy given by equation (3.14) may be written as 

U/(mn^R^ - I / (D^(p^ - e’)^ + + Dj. dx + | / EA(l-x*^^)^ dx 


+ i / (1 - x'*’^) {(EI^(p^ - e’’) D^pJ + °£;P;lx‘^^ - 2EAe^p^x‘^}dx 


+ 1 / EB^(p^ - - eb2Pj.(p| - e'^)x‘^^ 


-EOB 3 p^ Pj, + B^p^)x'^^ + - 2EC*p^(p^ - 0 ")} dx 

+ / G(J p + J p ) (p. - e')p.x‘''^ dx + U*/(tnfi^R^) (3.17) 

Q n n ^ s s s 

* 

The term U , which includes most of the contribution due to warping, is given in the 
Appendix . 

The quantities D^, and defined by equations (3.16a,b) are, respectively, 

2 4 

the torsional and flexural stiffnesses of the blade, normalized by mQ R , while 
— defined by the first of equations (3.16c) — Is the polar moment of inertia 

fi 4 

for the cross section (referred to its shear center C^) normalized by R ; A is the 

2 

blade's cross section area normalized by R , and e^ is the cross section area 
centroid of fset from the shear center, normalized by R. It is worth noticing 
that the curvature components p^x , p^x and p^x appear directly in the expression 
for the strain energy U. This Is a direct consequence of defining U via the strain 
components formulated In terms of the increments dr, dn and d(; for the deformed 
configuration. 
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CHAPTER IV - HAMILTON'S PRINCIPLE AMD THE DIFFERENTIAL 


EQUATIONS OF MOTION. AMP BOUNDARY CONDITI ONS 


In chls report the extended form of Hamilton’s principle (77,78] Is used 
to obtain the differential equations of motion for the blade, and their boundary 
conditions. This principle may be expressed as 


/ 


2 


T 


1 


(6(T) - 6(U) + 5W]dT 


A(nd2^R^)6l • 0 


(4.1) 


where $W — which is not necessarily the variation of a function W, and hence the 
Che distinction from Che notation used for iS(U) for example — designates the 
virtual work done by the external and damping forces. Here, 6W is expressed In 
terms of Che virtual displacements cf the non-dimensional generalized coordinates 
u,v,w and 6 as, 


6W 


{ 6W„ 


+ 


0 


1 

/ (Q 6u + 0 5v + 0^ 6w + Qq 60^)dx 

O ^ ^ X 


(4.2) 


The quantities Q , Q and Q are the distributed forces associated with the virtual 

u V w 

2 

displacements R6u, R6v and Rdw, respectively, normalized by mfi R, while is the 

2 2 

distributed moment associated with the virtual rotation 60 , normalized by mfi R . 

X 

The term 6W is included in equation 4.2 to account for the cases where the 

o 

virtual work done by the distributed forces acting on the system is also dependent 
on the variation of the spatial derivative of any of the variables u,v,w and 9^. 

In such cases, an integration by parts in the expression for oW immediately yields 
a 6Wg term, as in equations (6.15a,b) in Ch/»pter VI. 
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From equacions (2.3) and (2.8) to (2.10) the following quantities» which 
appear in the expressions for the variations of the kinetic and strain energies. 


are obtained 

• a • 

- ((9^ + cB)c6y - s0 s9y 00^)60^ - sB cO^ s0^ 60^ 60^ + 80^ 66^ 

• • • 

- t(6jj c0)s9y + tj^jSBlsOj^ 60^ - t22*0 00^ * ”*■ *23 

6u)^ • " ‘^ri tjj^sB]c0j^ 60y - t22s6 69^ + s0j^ 60^ tjj 60^ 


(4.3c) 

6Pc - e' c0 69 + 60' + 90 6e' (4.3d) 

C 2 y y X y z 

6p„ - p^ 60 - 0* 90 80, 60 - C0, 60’ + C0 90, 69' (4.3e) 

n Cxtyly ly ylr 

<5p^ - - p 60 - 9 ' 90 c0, 60 + 90, 60' + c0 c0, 60' (4.3f) 

C *^nxzylylyylz 

*'21 ■ '31 • '23 <=«. - '22 *®2 

*^2 ■ '32 ■ '23 ’*2 *®y '21 *®2 

6t-« • t,, 60 - 90, 90 60 (4.31) 

23 33 X 1 y y 

t 6x » - [(1 + u')6u* + v'6v* + w'6w']/[(l + u'^) + v'^ + w'^]^ (4.3j) 


( 4 . 

(4. 


The virtual variations 60 and 50 may be found directly from equations (2.11a,b) as 



69 . {. 


- w' 


ToTo^TTv^ 


l v^ 


* w 


• Oey/3u')6u’ + O0y/3v')6v' •► O0y/3w')6w' 

60 . [.v'* 6u' -► (1 + u’)«v*l/U u’)^ ■♦• v’*l 

z 

- (30 /3u’)6u* (30 / 3 v*) 6 v* 

z z 

From equations (4.1) to (4.4b) » (3.6) and (3.17) the following differential 
equations of motion are obtained after performing a few Integrations by parts in 
equation (4.1) 

{(Tg ■f Ug )30 /3u' (Tg + Ug )30^/3u' T^(l u') + h^^}’ ^ (x.X) 

y ' * t 

• ('i ’<■ ^ ®^22^ ^ + u -*• S - Qy 

{(Tg + Og )3e /3v' + (Tg + Ug )36 / 3v ' + T v' + h }' ^ c'(x,T) 

y ” * Z * ® '' 

-(v + eCjj)** 2(u + «t 22 )* cS • 2(w «t 23 )* »6 "(v 4- «t 22 ) ~ 0^ 

{(Tg + Ug )30 /3w' T w’ ^ h }’ L g;(x,t) 

y y y 

■(w + 6 X 23 )’* + 2 (v + et- 2 )* sS +(x + u + et 2 ,)s 6 c 6 -(w + et23)3^6 • 



^0 ‘ ^0 ‘'9 ^ % "0 

X X X X 


T ■ -X 
e 


‘ro^(p^ - e')^ + (3 - 2x*h (O^pJ + Dj;P;)/2 + ei5<l - 2”**) «>C - 


“Vc“ 


- 3x'^^)/ (2x'*’) + E(p^ - 0’2 )x‘*'[Bj(P? - 0 ^)x''' - SBjP^l/Z-EAd - x 


- BECBB^p^ Pj, ■*• B^p^)x"^/2 + 3G(J^P^ + J^P^) (P^ - 9’)P^x^ 




■*“ 90 c6y C0^] 


“ 

' 1^ ~^2v - ws6 +(x + u)c0^c^c9^c9y 4-1^’ + 2u cB- 2w s0 - ^s9^c0^ 

" 

^ C” ^ ^2v - w S0 4- (x + u) c6^s^s6y 's6j 


T - 30y + (Jr,<*)^s0^ + J^uj^c0j)c0y]* fJ^^cOySB^ + 


- j aj^(s9.c9 + c0,s9 s9 )]s9 - e^ 

s s 1 z I y 2 


[y 4- 2u c6 - 2w s6 - 


90 4- (x 4- u)c0] 


(A. 6c) 


30 




(Jp - " »9 ♦(« * “>'s) '^'31 


[* + Ji c6 - 2i la - ♦ [;• + ^2v - w gfl + (x "f u)csjs^cd^ 


6d) 


Ug -^ D ^( p ^ - e ’) ■*■ P^[(l - %^ hzi ^/2 ^ EB ^( p | - e ’^) x '^^/2 - EB ^ P ^ x ^] + ( EcJp^)Vx 


•♦■2 


* G ( JnP , H * J;P 


J ( 2 p , - e ’) x *^ x ^^ 0 ^ c 6 y - i (3 - x '^^) ( 0 ^ P ^ s 9 ^ D , p , ce ^) x '^^ e ^ s 9 y 


[it 


* j IIeAU - «*■ + EBj(p^ - e'^)/’ + lECBjP; + BjP*) 




t 


■>•3 


+ MEBjP ^ P ^ X •► 


EC *( p ^ - 9 ) Ja 0 ^ - G ( J ^ 30 ^ •*• Jj . c 9 j ^) p^(Pr - 9 ’) x '*'^ e ^ 


s9 
2 y 


^ |(3-x'*’^) ^ EB2(p| - 9 ^)x‘*‘ s6^ -t- 3E(B3P^ + B^pJ) 


x ^ sd^x 


H -2 


•f j EA(1 - *®l " 3 EBj PpP^x'^’’ c9j^ - ECJ(p^ - 9 )c9j^ 


+3 




♦ G ( J ^ c 9 j ^ - • 0 j ^) p ^( p ^ - 9 *) x "**^ 


T 


(4.6e) 
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•f f 1^(1 - + I Bj (p^ - e' •■ EBjP^P^x*^ + (EcJp^)V 


\*/ +2 


- |[d^(P 5 - e') + I 1^(1 - x“^^)p^ + I Bj(p^ - e' ^)p^x‘*’* •■ EBjp^p^x*^ + (EcJp^)Vx^ 

+ G(J^P^ + J^P^) (2p^ - 9')x‘^Jx'*'^ s9y + j(3-x‘^^’) (Dj^p^ «e^ f Dj.p^c9j)x^^c0y 

- i j^(l-x*^)e^x‘^ + EBjCp^ - e'^)x'*’^ ^BECBjp^ 

- [ 3 EB 3 -*■ EC*(p’ - 0")js03 C0y + G(J^ 80 ^ + Jj. c9^)p^(p^ - e')x‘^^ 0 y^ 


- P^®3 


A - H 


<4.6f) 


|jD^(p^ - 0') + f (1 - + f Bi(P^ ■ 9'^)P^’‘‘*’^- EB2 P^P^x"^+ (EC*p^)7x' 

♦ «->n% * '<^K> "‘’5 ■ * j”'**'* ‘V 


+ ifEA(l-/^)ey + E»j(p| - e'^)x-^^ + SECBjP^ + 


r +3 

3 EB 3 pj^pj.x + 


EC*(p^ - ®"j]^ 


(4.6g) 


The second term, EI_(1 -- x ) P^/2, that appears in the expression for 


Ug is a tension-torsion coupling term that reflects the increase in the effec- 

X 

tive torsional stiffness of the blade due to axial tension, e^S* [45,75]. This 
is discussed further in Chapter VII. 
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4 4 

Th« c«n& T^/x , glv«n by «quaCion (4.6«)i i« eh« eocffleitnc of -x 4x th«C 
appoar In eha variation of (U*U )/(on R ) givan by a<(uAClon (3.17). Tha caroa 
hy* appaar In tha lafc hand alda of aouaelona (4.5 a«d) danota 

tha amall contributiona dua to (T*- U*)/(Bfl^R^) . Tha taroa In aquation (4.1) that wara 
Integratad by parta yield tha following boundary condition aquation 

E C 5u + G 5v C 6w H 50 + E(C,0” - C?p„)60’ -EcV(a0 )«er + H f (30 /3u )«u' 

U V V XX IX inX lrv 7 ^ * ^ 

+ (30^/3v’)5v'J - Hy[(3ey/3u')5u' (39y/3v')5v* (30y/3w’)5w^- 5W^ -0 

X «0 

(4.7) 

where, 

®x * i i - 9'^)0^x‘^^ - EB^Pj-P^x*^ 

-fr ”*■ Ppi * 

and 

Hy -[y (3 - x**"^) (D^p^ c9^ - D^p^ 80^) + I EB2 (p| - 0'^)x*80^ + | E(B3P^ + B^pj)x'^s0^ 
- 3EB3P^Pj, x'*’ c9j^ + G(J^ c0j^ - 90^)p^ (p^ - 0')x’^x"*’^ 

+ I EA(1 - 90^ - EC*[(p^ - 9'’)c9^ - p^e; c9y) (4.8b) 
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CHAPTER V - EQUATIONS OF MOTION WITH O(e^) NON-LINEARITIES 

Equatio'As (A.Sa>d) arc a act of four non-llncar coupled partial differencial 
equations satisfying the boundary conditions extracted from equation (4.6) • In 
general, they admit non-r^iero particular solutions which are either an equlllbrlim 
state a ■ constant ^ (a • u,v,w,6^) for the caae of hover, or an steady state 
solution for the case of forward flight. Of special Interest Is Che determination 
of these particular solutions, and the analysis of the stability of the perturbed 
motion about them. Due to Che complexity of equations (4.Sa**d), they will now bs 
restricted to "moderately large" deflections. Toward this end, Che non-llnearlcles 
In chose equations are expanded In Taylor series and the results truncated to a 
certain degree in order Co obtain a set of approximate differential equations that 
are more amenable to analysis. It is well known that Che resulting equations, with 
polynomial non-llnearltles, still retain relevant information about the motion. 

Here the expansions will be carried to include terms up to third degree. 

To expand equations (4.5a-d) into third degree polynomial non-linearities a small 

parameter e, of the order of the bending deflections, Is introduced and an order of 

magnitude e"(n ■ 0,1,2,...) is assigned to Che variables and parameters of the system, 

e.g. [32,42,52,62]. For example, the blade's geometric pitch angle is taken to be 

of 0(1), but its space and time varying components [see equation (2.1)] are assumed 

to be small so chat 0 ^(x) •• 0(e), 6, • 0(e) and 9, ■ 0(e). The assumed order of 

pc IC 18 

* 

magnitude of Che quancicles that appear in equations (4.3a«d) are as follows 
Order £ : 


X, 0^ 





^ The order of magnitude of the quantities indicated after a semi-colon follow as a 
consequence of the assumption for the order of magnitude of the appropriate 
0(1), 0(g) and O(g^) quantities indicated before a semi-colon. 
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Ord€r e; 


’ ®x’ ®pt’ ®lc* *•’ *’ *A’ ®®2’ “ 3 * ®®4* 


Order C : 


'1^1 Jji^* 2Bj^, ECj 


In addlcloni EA • 0(e 

From equations (2.11a,b) and (4.4a,b) the 

for 0 , 0 and their partial derivatives with 
y z 

0 - w'(l - u’ - v'^/2) - w'^/3 0(£^) 

y 

6 - v’ (1 ~ u') - v’^/3 + O(e^) 

2 

00 /9u' - - w’ + O(e^) 

y 

99 /9v* ■ - v'w' -f O(e^) 

y 

90 /9w' - 1 - u' - v’^/2 ~ w'^ O(e^) 

y 

90^/9u’ - - v' + O(e^) 

90 /9v' - 1 - u' - v’^ + O(e^) 
z 

Also, from equations (2.8), (4.6a) and (2.1) 


following expansions are obtained 
respect to the elastic displacements, 

(5.1a) 

(5.1b) 

(5.1c) 

(5. Id) 

(5.1e) 

(5. If) 

(5.1g) 

it follows that 
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■ 1 - u’ - (v*^ >f w’^)/2 O(e^) 


(5.1h) 


I A (1 - K - EA{«» ^ I’ + f - 2lu* ^ f ^ f 1 } f 0(^") (5.11) 


.2 ,2 


.2 ...2 2 


• EA |u' 


r '2 „.2 u '2 

2 ”*■ 2 ■*■ 2 


f ♦ f " f ^ - W * 7 ^ - •jrv”a - u' «v’^ 


[[■ 


2 2 
y ) - V* -f ^w”(l - u' - ^\” W*^)“ W* (u" ■♦• v’ v")J*0j^ 


„.2 

- 6(u' +1 I ) (V 


ll\ 

" C0 -f w" 80 )> V - UvO' 

C cjj 

- ifD„(v" 80 - w" c6 0^(v’* C0 + w” 80 )MXo(e^) 

* L n c c ^ c ^ J \ 


|2 ^ i 10’ (0* + 20') 

,x 2 C x' X 


(5.11) 


Afccr axpanding che remaining terms In equations (4.6 a-g),‘«^nd noticing that 

\ It 

the terms h » h , h and hg In equations (4.5 S'-'d) due to the T an^ U 

* \ 

concrlbuclons to the kinetic and atraln energies are of higher order — \as can be 

readily verified from equatlone A.l, A. 2 and A.? In the Appendix'— the foll^^lng 


differential equatlone with polynomial non-llnearltles are obtained 


•f D_ (v" c6 + w" 80 )c9 - EAe.(u' + x 

s C C C A 4 


\ 


\ 

w'2 


O’ -K'rD„(v" 80 - w" C0 )80 

u 1 Ln ® c c 

r v*^ w'^ “1 * 

- w'l D_(v" 80 - w" c6 )c0 - D (v C0 + w" 80 )s0 EAe.(u' + x ^ )s0 1 

L,n c cct C cc A 2 ^ 

T 

I • • • • 

u )r • u - e( V 

- 2(v - e§, 80 )c6 + (w + e 80 )a0 c0 « | x + u - e(v' c0 + w' s0_,)X^6 

1 C c I C Cl 


- xe(v^ c0 w* sd )c p + T (1 + 
c c e 


c9 -f W 80 )* 
c c 


« f OCe-") 


(5.2a) 


2 2 

* Although slnS and cos0 are related as cos 0 -f sin 0 
COS0 Is left In the equations as c0 for simplicity. 


1, the 0(1) quantity 
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V + e(c0 



"(w* + 2v s$ 4- X s3 c8)c9 


J^(0j + w' c3) 4- (Jj.- J^) ^j:2v' 4-ce)89j^ c6^- w c 2 ejc 0 - •jCv 80^ - ¥ c9j^) (2v 

4- ( V* - 2w 80 - V 4- 2u c0>80 l.U;.-. 
-e;;(2u'4- v'^4-w’^)j4- El^ej (u’4 |'^4- l'^) - EB^ej (v” C0^4w”80^) j 
4.(v” 90^ - w" c9^) - D^) (v" c0^4- w" 80^) 4- EAe^(u’ 4- |'^4- |'^)J 


) (2v 4- X C0)c0 


+ (EC a")' - [EC,(v*» 90, - w*' C0,)3" - EC, 0”(v”c0 - w” 90 ) 

ii 1 IX C C 

- 1^(9;^ 4- ep 90^- w" C0^) 4 GJ^(v" C0^4 w” S0^)j| - Qq 4 0 


(e^) 


(5. 2d) 


In equaclona (5.1J) and 5.2 b-d) the expansions for sin0j^ and cos0j^, 
truncated to the appropriate order, should be used. With 0j^ expressed as 

epx.l) = 0^4 0 ^(x,t)4 6 (x) 4^ (6^^ ®c 

1-1 


(5.3) 


these expansions are 


90^ - S9^ 4 9^^ C0^ - i 0^^ 30^ - i 0^^ c0^ 4 0(eS 


- 0,, 90 - i 0^, C0 4 O(e^) 

11 C 


(5.4a) 


c9, ■ c0 


(5.4b) 


Equations (5.2 a-d) may be reduced to a set of three Integro partial 
differential equations In v, w and 6 , with cubic non-linearltles In these 
variables. For a rotor blade with the end x ■ 1 free, equation (5.2 a) may 
be Integrated with G^(x “ I, t-t) • 0 to yield an expression for u’ In terms 

of the remaining variables. To obtain such expression, the Integrated form 
of equation (3.2 a) Is first written as 


[D (v" s0, - w” c6 )]’ (w* C0 - v’ 80 ) x(v' c0 + w* 80 )e c^8 

ncc cc cc 


I 


+ (v' c9^ + w' 90^) j EAe^(u' ^ f ) “ D,(v" c6^ + w’ s0^) 


80 + (x -f u)c0 


- ej(v’ - v' c^0)c0^ + ( w ' - 20 j^c6 - s0 c 0 - w’ c^6)s0^J^ 


dx 


•d-u' ) / Q,, dx - •^’(l-x^)c^0 + O(e^) 


(5.5) 


Equating the right hand side of equations (5.1j) and (5.5) to each other, and 

multiplying the resulting equation by (1 6(u^ + v’^/2 w’^/2) - u ’/2 + (v’^ -f w’^)/4], 

2 2 2 

the following expression for EA(u’ v' /2 -f w’ /2), which is needed to 0(e ) in 
equations (5.2 b,c), results when Q • o(e) 
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2 2 

EA(u* Y I ^ ®®c ' *®c^ 


-rOpCv" c6 >♦■ w" s8 )V(v* C0 + w' 80 ) + X8(v’ C0 ■♦■ w' 80 )c^0 
Us C CJC c c c 

^ Vic^ ’ i 2 [v^'' «®c - ^ «®c '* *®c^] 

/^ii - ^2^ - w s6 + (x 'fu)e^eB - “ v* c^0)c0^ 

+ (w’ - 2§j^ c6 - s6 c6 ‘ w' c^0)80^^dx - ^ c0^ 


. f3 , ^ 1 / ,2 ,2 

[2 ^ 

^w" 90^ - V* w’ w” C0^ - V" C0j^ - v" 80J^ 


. 


- EA(v'^ •♦• w'^)^ - . / 4x + (EAe^)’ (v' c0^. + w' s9g)(u' ^ >♦■ 

+ Y (1 - x^) Q(u' + Y -f w’^)J c^S + o(c^) (5. 


6) 


The remaining u-terms in equations (5.2 b,c) are needed to O(E^) at moot, while 
2 

only an 0(£ ) expansion for u is needed to eliminate u and its derivatives in the 
right hand side of equation (5.6). To O(e^), equation (5.6) may be solved for u’ as 

u’ - e^(v" c6^ w" s9j^) " J (v' ^ + w’ ^) + -gj- [j (l-x^)c^S - f (2v c3 + Qy)dx] + O(k^) 

(5.7) 

Or, with the boundary condition u(x • 0, t ■ t) • 0, and making use of equations 
(5,4a,b), 
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U • c 6^ + w" 80 ) - 80^ - w” C0^) - (v’^ -f w'^)/2)dx 


X X 

+ / ^ (j(l ~ x^)c^3 - / (2v ce Q^)dxldx >*• O(e^) 


(5.8) 


When equations (5.7) and (5.8) are used In the right hand side of equation (3.6), 

2 2 

the expression for EA(u' + v' /2 + w' /2) becomes 


2 2 

EA(u* + ^’ 80 - w” c0 )TV C0^ - V* 80^) 

2 2Ln c cj c c 


- rO-(v” c9 4- w*' s 6 )T(v’ c0 + w’ 30 ) + xe(v* c0 + w' 06 )c^S ^ 0^.0 

LC c cjc- c c c ^ 

- i Elc0’ (0* 20') I- h)..(v” 80 - w" c0 )^ + D^(v" C0 + w" 80 )^1 

I H K c C 5 C ^ ml 

- / |^(2v - w 80 + X C0)C0 + (e-e^) J(v' -v' c^6)c0^ + (w* - w' c^S)s0 J 


,2 

X 


e(20^ + s0)c$ s0 




2 fl ,2 

c 0 - V 


- w'^lc^0 dx dx i(l -x^) (v'^ + w'^)c^0 " (v' ^ + w'^)^ 

J 


+ -|(3 c^0 - 1) (1 -K^)e^ (v" c0^ + w" 90^) + (1 - x^)^ cS 


+ EAe, 


V 

. .2 

“1 , ,2 

• 

v" 

< 

N 

+ 

n>|€ 

) ■♦■ v' w' w’Mc0^ + (w*'^ + I )w" 

S0 

c 

J 


/ Q„ dx 


+ EAe. (v" c0, + w" 80.) - -r EAe?(v" c0 + w" s0 )“■ + o(e^) 
A1 12A c c 
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(5.9) 


The integro partial differential equations In v, w and 9 ^ with non-llnearltlea 
to 0(c ) are obtained when u and Its derivatives In equations (S.2b,d) are re- 
placed by the expressions given by equations (5.7), (5.8) and (5.9). The re- 
sqltlng equations are written here as equations (5.10) to (5.12) for the special 
simpler case of constant stiffnesses and e * e^ ■ 0. 


w*^) 


•i T[% - !" - “>n \ 

* Ir [“ - I <1 - - <®n \ [l 

+ [v - I (1 - x^)v"J*j -|^w‘ *► v” w’) + EcJ(v”’ s 9 ^ - w” ' c 9 ^)j 

- I '*’ 29 ') + 3 B 3 (v" s 9 ^ - w" c 9 ^)^ + 3 B^(v" c 0 ^ w" 89 ^)^Jc 9 g 


- 3 EB,(v" s 9 - w" c 0 ) (v" c 9 w" s 9 )s 9 - EC, ( 9 * + v” w')' s 9 , 

3c c c cclx I 


+ G(J 


■f X c0)c0 


„ s9 *«■ J, cSJQl (9' + 9')1 - v'f / f(2v - w s0 

ijC ^CXX I fj^W 

|||- (x - J “♦• / / V* + w' w*)‘ - -J (v’^ + w’^)c^^dx dx 

[^v’2 + w’2 - (1 - x^) j + s9^ - w c9^)’ s9^ 

+ J^(V c9^ + V s9^)' c9^ + ^^^1 ■*’ *®c ^®cl 


^ s . 1 ^ .I 

- V - 2 w 30 - 2(c0) / |v' v' +w’w» + — / (2 V c$ + Q^)d^dx - 4 . O(e^) 


n 2 . c ^0 

- (1 - X ) -j- 


s • 

V 


(5.10) 
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G’ - 
w 




'“n - "«) I *Wi - <">0 ' «1 * • « 


*f (D^ - w' w 


tf# 


^ [l ■*“ *' 1^ ‘ 2 " *^)v'’J*j .8^ c9^ 

■ % *^®c ^ °c [w - I (1 - x 2 )w"]’ 

+ w'(w’ w")’|4. v” [o^(9i •*■ v" w*) •► EcJ(v'” ge^ - w”’ C0^)J 
+ ^ 82 ®; ( 9 ; + 29 ’) + BBjCv” i 9 ^ - w" c 9 ^)^ + 3 B^(v” c 0 ^ + w" » 9 g)^* 9 c 


- 3EB-(v" 80 - w" c0 ) (v" c9 + w" s9 )c0 - EC, (9' + v" w’)' c9, 

JC C C CCiX 1 


+ G(J^ c0 


c'-’c *v*; < 9 ;*«')] - »'f/’'[( 2 ;-».« 


+ X c8)cS 


■*■ flf “ f ) 0*^’' ■*■ / / Rv’ v' + w' w') - (v’^ + w'^)c^8*jdx dx 

J 1 0 I- J 

- (1 - n.h ^ + w'^ - (1 - (V 39 ^ - i c6^)' c9^ 

■f J^(v c9^ w s9^)’ 80^ ‘ (J^ - c20^ c6 - j^(9j + s8 + w' c8)c8 




+ 0„. s^9 + c^9 ) (w' c8 + s8)c8 V - w + (2v - w 

" *i c c c 1 


s8)s8 


J 


[x - -j / (v'^ + w*^)dx + (x - c8 - O(E^) 


(5.11) 
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w’ c8)* (j^ - j^) Qv’ + 1^) 826^ - w’ c29gJcS 

‘|°5[®x ‘ - x^> + EI^9[(1 - x^) - EB^e’^Cv'' c6^ + w” 80^)^ 

+ (D^ - Dj.) (v” 80J - w" c9^) (v" c0^ -*• w" 80j^) + EC^e”" - EC*J(v” 80^ - w" c0j^)" 

9'’(v" c9 + w” 80 )1 - 1 ( 9 ’ •► 0 ’) rSJ«(v" 80 - w” c0 ) 

X c ^l/^^Ln c c 

GJj.(v” c9^ + w" •0c)jj 


Qa + 0(0 


(5.12) 
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CHAPTER VI - GENERALIZED AEROD\TJAMIC FORCES EXPANDED TO O(e^) NON-LINEARITIES 

The velocity of a point on the blade's elastic axle, relative to the hub 
center 0, may be obtained directly from equation (3.2b) In Chapter III as 

v^ • QR{(u’- V c0)i (v + (x + u)cB - w 80] y + (w + v s0)z } (6.1) 

A 

The rotation of the blade, with angular veloclfy HZ (see Fig. 2. 2 in Chapter II), 

A 

Induces a "downwaah" air Inflow with velocity - HRv^^Z • - nRVj^(ft s8 + z c0). It 

should be noted that In the notation used here v^ denotes the magnitude of the 

Induced Inflow, normalized by QR. If denotes the magnitude of the constant 

forward velocity of the hub center, 0, normalized by f2R, the velocity v . of 

6 / a 

the blade's elastic axis relative to the air Is written In terms of the advance 
ratio p A cos ■ 0(1) and Inflow ratio A A sin "*■ * 0(e) as (6,9,28] 

A A A 

V / ■ v + flRV,(Xca + Zsa ) + DRv.Z 

e/a e f r r 1 

- v^ (2R{p((ii c0 - z s6)cT - p sx] + X (x 80 z c0) } (6.2) 

Making use of equations (2.2) and (2.3) In Chapter II, equation (6.2) may 
be expressed as 

AAA 

”e/a " ^T^ (6.3.1) 
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where 


m m 



\ 


u - V c0 + u c$ cT + X s0 

s 

- [T] 

V (x <f u) cB - w s0 - M sr 



w >f V s6 - y s0 cT X c6 



— ^ 


The aerodynamic forces and moments will be modelled using Greenberg's 
extension of Theodorsen's theory in which only Che and Up compone*'ts of the 
velocity are assumed to effect the aerodynamic loading [17-19,34,46,52,56]. 

According to this theory, the lift L and aerodynamic moment M per unit length 

^ A 

along the blade, expressed In the (5,n,^) unit vector triad, may be written for 
quasi-steady aerodynamics as 

L - (mH^R) (L^(; ca + n so) +• (6-4) 

- (mn^R^)M^5 (6.5) 

where the non-dimensional quantities and are given by [56] 


* A more general expression for the generalized aerodynamic forces In hover may be 
obtained by allowing the aerodynamic coefficients to be a function of angle of 
attack [82,83]. As shown In [83], this procedure yields a refined expression for 
calculating the steady Induced inflow velocity in hover. 


h 


47 


1 


Sc - 


24 S ^ 


(6.6b) 


S • * 6 f *A ""S ^ ■ S>S "■ !<il - "*A ^ 2 Ka)«^1 


As illustrated in Fig. 6.1 the circulatory lift component mil RL^ acts in the 
direction normal to the velocity ilRU , where 


U - (uj + 


(6.6c) 


(6.7) 


The angle a is the angle of attack for the airfoil, while c and are, 
respectively, the airfoil chord and its aerodynamic center offset from the 
shear center C^, both normalized by R, These quantities are shown in Fig. 
6.1. The non-dimensional parameter Y that appears in equations (6.6 a-c) is 
the Lock number for the flow, defined as 


Y • bxp.R c/m 


( 6 . 8 ) 


and p is the air density. The distributed aerodynamic drag acting on the blade 

is also shown in Figure 6*1 and is written in terms of the airfoil profile drag 

coefficient c as 
(?U 

®a “ y (C sa - n ca) - (mfl^R) J U^(C so - n ca) 


A (mJl R) D (C sa - n ca) 


(6.9) 


As can be seen from Ftg. 6.1, sin a and cos a may be expressed as 


sa • - Up/U ; ca • U^/U 


( 6 . 10 ) 
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Cooiblnlng equation* (6.4), (6.9) and (6.10) the r**ultant, of the 
distributed aerodynamic force* acting on the blade 1* deteimlned a* 


P^/(mn^R) - (L^, ea - D co)n + (Lj,j. L^. ca + 0 •«); 


S Fr, n ♦ 


< 6 . 11 ) 


where 


f„ ■ I <f - ’'A>"p “5 - ^ 


(6.12a) 


■ 6 I- “p“t * <! - -A* Vs - f ”p * 


«a>“5 - # "“p> 


(6.12b) 


The virtual work done by during a virtual displacement R(x 6u + y 5v + 2 5w) 
Is then 


(fiw)^ . = - (mfiV) / (F n F, C) • (x 6u + y 6v + 2 

'due to F^ 0 ^ ^ 


cSw)dx 


- (ni7^R^) (L 5u + L 6v + L 6w)dx 

' M V W 


(6.13) 


where, according to equations (2.2) and (2.3) In Chapter II, 


“ t,, F + t-, F 
u 21 n Ji ^ 


(6.14a) 


^v “ ^22 *^32 


(6.14b) 


L ■ t,- F„ + t-. F 
w 23 n 33 C 


(6.14c) 
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1 


Wich eh« virtual rocatlon about tha ^ “ axis of tha airfoil obcaload dlractly 

A 

from althar aquation (2.9) or (2.10) at [39 + tha virtual work 

K y z 

aasoclatad with tha aarodynamlc oomant la 


m 


dua to 


- n (rfjV) / Mj«0 (06 /3w’)6tt’ -► oe./6v’)dv’]a0„}dx 
LAX z z y 


(mO 


3 r I* I 

«H, + / 

L "» 9 


{M^60 - (M.06 /3u’)a0,l’6u-(M.O0 /3v')a6^]'3v)dx 

Ax A s y A S y 


] 


(6.15a) 


whara» upon aubatltutlon of a0 by v* x*, 

7 


6Wa - M. (09 /3u')3u + 09 /3v')3v]w»ir 

IS A Z Z 


(6.15b) 


'i 

.1 


The ■'SWg term given by equation (6.15b) Is Included In the boundary condition 
equation for the system [equation (4.7) In Chapter IV]. 

Equations (6.13) and (6.15a) yield the following non*dlmenslonal general- 
ized forces due to the aerodynamic loading 


^*^u\ero 

L 

u 

- (M^O0^/3u’)w'x'^]’ 

(6.16a) 

^'^^aaro 


- (M^O0^/3v’)w'x'*’]' 

(6.16b) 

aero 

"w 


(6.16c) 

(Qfl " 

u aero 

«A 


(6.16d) 


5l 


If there are no other forces acting on the blade, such as dani|plng tor v«xampla. 

To expand the generalized aerodynamic forces Into polynomial non^llnearltlas, 

the angle of attack, a, of the airfoil Is assumed to be small so that Up “ 0(e). 

With ■ 0(1), equations (6.12 a,b), (6.14 a-c) and (2.3) Indicate chat the 

3 2 

expansion for Up Is needed to 0(e ) and Che expansion for to 0(e ). We then 
express Up and U^ as 


Up - Upj -h Up2 + Up3 0(e") 


(6.17a) 


“t • “to * “ll * -12 ♦ > 


(6.17b) 


where and denote the order c. (1 * 0,1,2...) term in the expansions 
for Up and U,^,, respectively. With 


y/6 - 0(1) 


c - 0(e) ; - 0(€) 


c^^/(2n) - 0(e ) 


and 0)^ ■ *f w’ c0 s(5 4- w’v’ ^ O(c^), as given by equations (2.9) and (5.1a,b), 
the following expansions are obtained from equations (6.12 a,b) for and 




(6.18a) 


^ri ^2 + ^ > 


(6.18b) 
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where 




( 6 . 18 c) 


n3 


ig 

(Y/6) f2Upj Up2 - (f - x^) Upj (5j + w’ cS sS) - - O(e^) (6.18d) 


‘^Pl ‘^TO “ 


( 6 . 18 e) 


C2 


(Y/6) {Up^ > Up 2 U^Q - (| - x^) (§^ + w' c8 4- .6 )U^q + f Upj) - O(E^) 

(6.18f) 


F - - (y/6) {Up^ U^2 ''P2 "n * "^P3 “to ’ -*A> + w' cS + 80)0^^ ^ w-VU^^) 


^ f “P2 - f<f - ^ ^ # “to “Pl^ ■ 


Making use of Che expansions for slnO^ and cosOj^ given by equations (5.4 a,b), 
the following expanded forms for the generalized aerodynamic forces are then 
obtained from equations (6.16 a-d) , (6.14 a~c) and (6.6 c) 


(Q ) 

u aero 

- (v' 30 - w' c9 ) F^, + 0(e'’) 

c c 

(6.19a) 

(Q ) 

V aero 

■ - »®= * «’'n2 - «11 ’’a>'*= - "u 



* '<^3 - «u ^2 - 

- "na - i * ”a - f 



+ <Y/6)x^[w'U.j.qUp^]' + O(e^) 

(6.19b) 
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"a * 'C2 '»c> 


2 

■*■ ' ®ii *’c2^*®c * ^®ii *"n2 ■ i ®ii *"f: 3 “ 2 


+ 0(eS 


(6.19c) 


^*^0 ^aero 

X 


6 'V'Wpi * "io“p 2 * "ti"pi> + < f - *' 'S * . 6 )U 

f< f ■ ^"pi’ * 


TO 
(6.l9d) 


To complete the formulation of the generalized aerodynamic forces, only the 
closed form expressions for the terms Up^ and (1 “ 0,1,2,...), as defined by 
equations (0.17 a,b), remain to be determined. Making use of equations (2.3), 
(5.1 a,b) and (5.3), the following approximate expressions are obtained for U.p 

and Up 


U.p = (v' c6j^ -«■ w' s6j^) (v - M ct)c 6 + (v u c0 - w S0)c0^ 

,2 

+ (x C0 - p ST) f(l - Y - v' w' s9j^] + (w +(v - M cT)s 6 + X cB)s9^ 

(6.20a) 

,2 

Up a (u + X 90 - V cB) (v' S0^ - w' c9j^) P(cT cB) fv' (1 - u' - ^ )s9j^ 

2 

- w'd - u' - v' - J )c0^1 - (u cB - W 9 B)s9j 

,2 

- (v - y ST 4- X cB) ((1 - Y + V sB c9j^ 

w'2 

4- (w - y cT sB 4- X cB) (1 - y )c6j^ (6.20b) 
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with 96^ and expressed as In equaclons (5.4 a,b), and with Up when 
u - " - « - 011 - 0 being an 0(e) quancley, che following expressions for 
U^^ and Up^ (1 ■ 0,1,...) are obtained from equaclons (6.20 a,b) 

- (x cS - U 9T)c0^ (6.21a) 

U-, ■ - (v' c0 + w' s0 )u cT cB + V c0 - (x c8 - U st) 0, , s0 * 

c c c 11 c 

+ (w - U CT s6 + X C0)s0^ 

(6.21b) 

U_„ (v' c9 + W* 90 )v cB + (v’ 90 - w’ C0 )0,, li CT cfl + (u cB = w sB)c0 

c c C C 11 c 

+ (v sB - V 0 )30 - (x cB “ li 8T) (“ C0 + v' w' S0 + •“ 0^, C0 ) 

11 c 2 c c 2 11 c 

+ (w - u CT sB + X cB)0,, c0 (6.21c) 

11 c 

Up, “ [ (X c0 - y CT sB)c0 - (x cB - y 8T)30 ] - (x c 6 - y st)6,, c0 

rx C C 11 C 

+ (v' 80 - w’ c0 )y CT c6 - V s9 + w c9 (6.21d) 

c c c c 

Up 2 - (w' c9 - v' 80 )v cB + 0,, (v' C0 + w' 80 )y CT cB - (u cB - w 9B)80 
^ ^ CCilCC c 

2 

- V 9 c9 + (x cB - y 9T) s0_, - v' w' c9 + ^ 0?, s0 ] + v sS c9 

C X C C Z 11 c c 

- (w - y CT sB X cB)0., s0^ (6.21e) 

11 c 
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U - - (u + A 30) (v' 30 - w' c9 ) - 9,, (v* c9 + w’ 30 )v c6 

C c ii c c 

2 

- U(CT c8) (| 9^^ <v' 30^ - w’ c9^) + v'(u’ + ^ )s9^ -w' (u' + v’^ -h 

- (u c6 - W S0)9j^j^ c9^ + ” v’w* c9 ] 

2 

•*• (X c6 - U 8T) [(| 9j^ + I’ )C9^ + V’ «• 30je-^^ - 9^^ V s8 s0^ 

- (w - y cT s3 + X c0) (9^^ + w’^)c0^ (6.21f) 



1 
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K>|C 


A case of particular interest, for which the equations of motion reduce 


to a simpler form, is that of an inextensional blade. For this case no EA 

terms appear in the differential equations of motion as the blade behaves as 

Infinitelly stiff to axial extension. The equations of motion for this case 

may be obtained by taking the limit EA ® in equations (5.2b - d) , with u 

3 

obtained from equation (5.8) in order to recover the 0(e ) inextensionality 

constraint [ 55 ]. It should be noted that the spatial independent variable for 

this case is now the arc^length r(0 £ r £ 1) along the blade’s neutral axis. 

A question that naturally arises is ’’what line along the blade’s span is 

the inextensional line.” It can be shown with the aid of equations (5.7) and 

(3.9b) that in the limit as EA the line that connects the area centroid 

of each cross section (point C (^ - 0, n ■ e*, ^ ■ 0) shown in Fig. 6.1) along 

a A 

the blade’s span is the inextensional line. According to equations (3.9b) and 
(2.7a ^ c) , if the line of area centroids is inextensional, the following 
condition has to be satisfied when warping is neglected, for simplicity, 


^ * [(1/x - f (I/>^ * ^12 ^ 


32 - 


+ Ki/x'" - 1 


or 


{ f(l + u') 


9 9 

+ u')“ + 






1 9 9 

} ‘ + (e.Pr)“ - 1 



It can now be readily verified that, in the limit EA equation (5*7) is 

3 

simply the 0(e ) expansion to equation (7.1b). 

An alternate way to obtain the differential equations of motion governing 
the flexural - flexural - torsional motion of an Inextensional blade via 
Hamilton’s principle is to adjoin the inextensionality constraint to the 
Lagrangian of the motion by a Lagrange multiplier as in (62, 68). The advantage 
of using the limiting process previously mentioned is that both the extensional 
and inextensional cases may be investigated via the same equations of motion. 

The differential equations of motion for an extensional rotor blade were 
derived in this report in terms of the three elastic displacements u(x, t) , 
v(x, t) and w(x, r) for any point on the blade’s elastic axis, and an Euler angle 
. The total angle of twist of the blade, obtained by 

first expressing the twisting curvature defined by equation (2.10) 

as (32, 33, 55, 61, 63] 

= x'^(0' + 6' + e' sP ) A x'^p^ (7.2) 

t L X a y - C 

and Chen integrating equation (7.2) over Che domain 0 £ x £ 1. Since the total 
angle of twist, , includes the non-elastic pitch angle P(x, r) , the elastic 
angle of twist of the blade, <J)(x, t) , is then 


(p(x, t) - / ('ll’ - 0’)dx + ({) (x=*0, t-t) - 0(x-O, t-t) 
0 


(x, r) + / 6 ' s0 dx = 0 + / v' ' w* dx + 0(e ) 

‘ 0 ^ ^ 0 


(7.3) 


Equation (7.3), relating 0^ to (j>, could be used in Che outset to derive the 
dif ferenti.al equations of motion for the blade in terms of u, v, w and tp, IThen 
this approach is taken, Che expression for i5p^ for example, given by equation 
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(4.3 d) , siraply becomes , In such case, the term D-(p- - 9’)5pr • 

s s s s 

obtained when the variation of the strain energy U given by equation (3.17) is 
taken, gives rise to a single term D^(|)’‘ in the 6 ^ equation. In contrast, If 
-0^ is used instead, the term 6p^ gives rise to the term in the 

60 ecjuaclon but also to the 0(e ) terms (D.iji'w')'' and (D-i^*v’*)' in the 6v and 

X S ^ 

6w equations respectively [equations (5.2 b, c)). However, the differential 

equations derived via 6u, 6v, 6w and 69 as done here and in [62, 63), and those 

2 

derived via 6u, 6v, 6w and 6({), ats done in (32) to 0(e ), are equivalent to each 
other since one set of equations may be converted to the other by the integral 
transformation given by equation (7.3). 

■f2 

As mentioned in Chapter IV, the term ET^(1 - x )Pf/2 that appears in 
equacion (4.6 g) is a tension - torsion coupling term that reflects the increase 
in the effective torsion stlffneas of the blade due to axial tension. To show 
this, the case of a non-rotating beam undergoing only static axial extension 
and torsion, as in [75], is nov/ briefly addressed. For simplicity, it will be 
assumed that no external torsional moment, , is applied to the beam and 

X 

that the beam’s cross section is symmetric about both the n and ^ axes. From 
the expression for the variation 6U of the strain energy, the following 

•k 

differential equation is obtained when the coefficient of 60^ is equated to sero 

D^_e' + I 1^(1 - x'^'^xe’ -*- 0') + I B. (0*^ + 20'e')(0' + 0')x‘^^ 

- (EC^9^’) Vx"^^ > 0 (7.Aa) 


Here D-, for example, denotes the torsional stiffness of the beam, as defined in 
Chapter III, but without carrying out the normalization process implied by equations 
(3.15 c) and (3.17). A similar observation applies to the other quantities in 
equations (7.4a,b) • 
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X 

With * 1/(1 + u^) - i - and u* ■ - (/ Q^dx)/(EA) ■ T^/ ^EA) au given by 
equations (j.7) and (5.5), equation (7.4a) may be approximated as 

(D. T I./A)9’ > T I. 0VA - (EC, 9'’)’ 4- O(c^) - 0 (7.4b) 

^e^xe^ lx 

The coefficient of 9^ in equation (7.4b) is the effective torsional stiffness 
of the beam [14, 75, 84]. For slender beams, the first two t nns in that 
equation are the dominant ones and equation (7.4b) reduces to equation (18) in 
[14]. A more accurate formulation for the torsion - tension coupling problem 
that cakes into account the non-orthogonality of the coordinate system (x, n* 0 
for the elastically undeformed pre-twisted beam v^;as recently presented in [75]. 
The work in [75] diw»ciosed that an extra term proportional to 0’ should be 
present in equation (7.4b) and that for values of G/E that are typical of 
isotropic materials the influence of this extra term is quite small and may be 
neglected. 
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APPENDIX - KINETIC AND STRAIN ENERGY CONTRIBUTIONS t! U* 

■ — - ■' ™ ^ ■ ■■■ " ' ' ' ■ ■ ■ ’ ^ 


The warping terms not shown explicitly In equation (3.16) are Included in 


T which is given as 


* 2 1 
T /(mfi^R*’) 


0 ^Tl^fl2‘^C* ‘ ^^ll^C^ ^,'^31 ‘ 


- [^12 ^^33 ^ 


32 ^31*^® ” tjjSS) 


/ S2 ^ Qi2''^ 

+ [(t^3sS - t^^c8)p^x^ - (t^ 2 Pc’‘'^>‘]^j' 


ti2P^x-^s8 + (t^3P^x+)^^ 


where 


mC^^ - R^ // ptC dn dC ; mC^2 “ 

A A 


'Fhe expression for the term appearing in equation (3.14) is 


^ R r 

“s-2 i 




(np^ - cp^)9iij/9n - 'i'P^T P^x'*’’^ + ^(np^ - ;p^)9t|;/3; + ii/p^ p^x'*"'^ 




+ III OK^/3n)'^ + (WdQ) 


2^(d\ii/dO 


fp.(p.x'^) 'x 


- e’e'T {* (n-91^/90 (p,- - 


'*’ - e'e'H I (n-9i|j/9c) (Pf - e') + f(np.-cp„)9ijj/3j;+ii^p. 


- 2ij>(9'l'/3n) |j)^(p^x'^) - 6 8 J j^(c + 3<i^/an)(p^-e ) + j^(np^-;p^)3^/3n-'i^pyc 

+ - 9'^ (3'»'/3n) (?<!'/ 3r,)J^dn dcdx 


The strain energy contribution U appearing in equation (3*17) is obtained 
directly from equation (3.14) as 


it it 

U « U +< 

s ^ 


i / //e (L^ + H 3 ) + (H ^3 + 2 L 2 H 3 ) 


(H^ + L^) + 4a^ ■'■ ^l”3 j 

|/ //E [a 

0 A L 


'dn dC dx 


4 2"^ 2 

(cp^ - npJ + 2(n + C) (i;Pr, - np.)^' (p^ 


- e ) 


dn 


+ 1 /eCj - e" I 

- /EcJp^j^^PjX^'' - 0"x*]x 


- (p' - e")Jdx 



2 q 

mil 


where , 


2 

4.9 4. 2 4-2 2 2 2 ’ -f 2 

- i - x^“ - 2 np^x^ + (cp^ - np^) X ^ +(n^ + ;; ) (p^ - e )x 


2C P„x 


,2 


1 I ~ 9 2 2+2 

[(p.x^) + (p; + ppp^x 

S ‘ I S *3 


„2 

- 0 ]x 


+2 


" +'’ + ' + 
» 0 X “ - (Pr>t ) X 


2 4*3 4“ ^ 4*2 

« (cp^ + np^)p^x - (CP^ - HP,) (p^x ) X 



(A. 2) 


dp dx 

(A. 3) 

(A. A) 
(A. 5) 
(A. 6) 
(A. 7) 


(A. 8) 
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